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Abstract 

Given an iR'^-valued supercritical branching Wiener process, let 
'il){A,T) be the number of particles in A C M"^ at time T, (T = 
0, 1, 2, . . .). We provide a complete asymptotic expansion of V'C^) ^) 
as T ^ oo, generalizing the work of X. Chen (|2j). 



1 Introduction 

Consider the following model in IR^ (with d > 1): 

(i) a particle starts from the origin in IR^ and executes a Wiener process 
W{t) e R'^, 

(ii) arriving at time t = 1 at the new location it dies, 

(iii) at death it is replaced by Y offspring where 

P{Y = i}=pe, (£ = 0,1,2,...) 

oo 

1 < ^ ipi = m < oo, 

oo 

< ^{i - mYpi = (j^ < oo, 

£=0 
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iv) each offspring, starting from where its ancestor dies, executes a Wiener 
process (from its starting point) and repeats the above given steps 
and so on. All Wiener processes and offspring-numbers are assumed 
independent of each other. 

Let 

J 1 if a; G IR^ is occupied by a particle at time t, 
I otherwise. 



A(x,t) 
We write 



which stands for the number of particles at time t located at A C IR^. In 
particular, %p{]R'^,t) is the total number of particles alive at time t. 

Since the branching is supercritical, it is well-known (Athreya and Ney 
d], p. 9) that 

(1.1) iVo := lim a.s., 

exists (and is finite), and that 'P{Nq > 0) > 0. 

The limit properties of ip{A,T), T ^ oo, were studied by Chen (|2]) who 
proved 

Theorem A. There exist random variables Ni and N2 (Ni being iR -valued) 
such that for any Borel set A C IR^ with ^4 ||a;|p(ia; < 00, we have, almost 
surely when T ^ 00, 

(2vrT)"/2^KAr) ^^r _ _L f ^No\\xf - 2N^ ■ x + N2)dx + o(T-i). 

J A 21 J A 

This result plays an important role in Revesz (0) in the study of the 
concentration of particles in the branching process. 

The goal of this paper is to provide a complete asymptotic expansion for 
'ip{A,T)/m'^ as T ^ 00. Let us first introduce some notation. 

If a = (ai, . . . , aa) G Z't and x = (xi, . . . , xj) G R^ we use the notation 



\a\ 



ai H had, a\ = ULi 2;" = nf=i x^' and 



(1.2) M„(A) = / x^'dx. 

J A 

If also (3 G Z'^ we will write (3 ^ a to mean that (3i < ai for all i, and if 

/? ^ a we set 

Here is the main result of the paper: 
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Theorem 1.1 There exist random variables {Na, a G Z'l) such that for any 
k > 1 and any bounded Borel set A C iR'^, when T oo, 



(1.4) (27rT)^/2^^^^^ 



E^P E ^ E g)(-i)""A«i)iv.o-, + o(r- 



a.s. 



n=0 ^ \a\=n /3d2a 

Remark 1.2 The random variables (A^^q, a G Z^) are described in the proof 
of Theorem ll.il They are limits of explicit martingales related to the branch- 
ing Wiener process. 

Although the distributions of the random variables (A^^, a G Z^) are not 
known, Theorem 11.11 can nevertheless be used to make predictions to any 
degree of accuracy. 

To see this, choose an integer k and disjoint sets {Aa C M'^, \a\ < k). 
Consider (1.4) for each Aa. Then we have a linear system of equations with 
the unknowns N2a-p- One can solve this system of equations if the corre- 
sponding determinant is not equal to 0. It is easy to see that we can choose 
the sets A^ such that the determinant is not for any T (T = 1, 2, . . .). Ob- 
serve the number of particles of a branching Wiener process which are located 
in the above given sets (A^, \a\ <k) at time Tq. Having these observations 
one can evaluate the actual values of the random variables {Na, \a\ < k) 
with an error term o{Tq^). Having these values one can use Theorem II. II to 
get the values of the process (27rT)'^/2^(/l, T)/m^ for any A C iR-^, T > Tq 
with an error term o(Tq~''). 

The proof of Theorem II. II is presented in Sectional In Section El we show 
that if the offspring distribution Y has p moments for some even integer 
p then the martingales described in Remark 1.2 converge to the random 
variables {N^, a e Z^) in Lp. 

2 The proof 

We start with a preliminary result concerning the transition kernel of the 
Wiener process. Let 

^^''^^)=(2;^^"P (-^)- 
Define the Hermite polynomials by 



(2.1) H^{x,t)=J2 



j!(^-2j)! 



3 



Lemma 2.1 For any < t < T and any x G M^, 



(2.2) f-W=(^|^?^^-(-.')- 



Proof. Let us recall the Hermite polynomials: 

[n/2j 

(2.3) = «! E ^^^(2.)"-^ . e 

so that 

(2.4) H„.{x,t) = {t/2r/^H^(^y xeR, t > 0. 



We use the following identity, see for example Lebedev ([3], p. 75): for 
any a > and y & M, 



2, ,2 11^ 



e 



;-i)v 



5o22"n!(l + a2)«+(i/2)^2n(y). 



Taking y = xj y2i G iR^ and a = yt/{T — t), and multiplying both sides by 
(27r(T - ^))-^/^ we readily get (jOI). 

If a = (ai, . . . , ad) G Z'l and x = (xi, . . . , Xd) G -R'' we use the notation 

d 

(2.5) i7„(x,t) = n^c..(x„t). 

1=1 



Lemma 2.2 For any < t < T and any {x, 

1 CO / 1 

^ n=0 |a|=7i 

and 

(2.7) pi.l(x - y) 

^(^l^,£i§„f;)<-''^---<-"- 
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Proof. Since for x = (xi, . . . , Xd) € 



(2.8) pi'\x) = l[pi'^' 



Xi 



i=l 



()2.fij) follows from ()2.2j) . To obtain ()2.7j) we use the fact that 



n ' 



(2.9) Hr,ix + y,t) = J2[ .\x--^H^{y,t). 

j=0 \J / 

For this we recall that (Lebedev 'Sj, p. 60) 

oo n 

(2.10) ~ 



[x] = e 



n=0 



so that 

(2.11) Y.^-lHn{x,t) = e'^-'''/'. 



oo 



n=0 



Then 



oo n 

(2.12) -Hn{x + y,t)= e^(-+^)-*^V2 

n=0 

oo A; A: oo j 

k=o j=oJ- 

and comparing powers of s"' proves ()2.9|) . <C> 

Now we turn to the study of the branching Wiener process. Clearly, for 
any T > 1 and A C R"^, 



E 



(V^(A,T)|^(r-i)) = m/ Ypi\y-^)Ky,T-i)dx, 



(as usual, J-'(t) denoting the cx-algebra induced by the branching process until 
time t). A simple argument by induction yields that for all < t < T, 

(2.13) E{^P{A,T)\m) = m^-' [ YptUv ' x)Xiy,t)dx. 

J A y 

It turns out that ip{A,T) is quite close to its conditional expectation, as is 
confirmed by the following results. 
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Fact 2.3 (Revesz |4j, (6.16)) Fix 7 G (0, 1) and let t = [T^J . Let A C 
he a hounded Borel set. Let e > 0. We have, almost surely for T 00, 

(2.14) ti^-l-J Yp?-tiy-^)Hy,t)dx = o(m-'/('+^^' 



Fact 2.4 (Revesz [4 , (6.11)) There exists a constant C = C{m,d) > such 
that for all 1 < t < T , 



(2.15) E y: {My, T) - E [A(y, T) I m]r < c 



^2T t 



(T-tY/2- 



Lemma 2.5 Let e > 0. Almost surely for all large t, we have X{y,t) = 
whenever \\y\\ > 

Proof. This follows from the usual estimate for the tail of the Wiener process, 
the Borel-Cantelli lemma, and (jl.lj) . <0> 

Lemma 2.6 Let a G Z^, and let 

(2.16) V^{t)=YHa{y,t)X{y,t). 

y 

Then, {-^Va{t), t > 0) is a martingale and 
exists and is finite almost surely. 

Proof. We start by proving the martingale property. Recall that '4){IR'^,t) 
stands for the total number of particles at time t. Thus, by numbering these 
particles and considering them all starting from time t = (many of them 
share common paths, at least partially), we can write J2y Ha{y,t)X{y,t) = 

J2fj^f'^^Ha(W^i)(t),t), where {W^'\ 2 > 1) is a sequence of iR-^-valued 
Wiener processes (they are not independent). Conditioning on jF(t — 1) 
and on iplM'^, t), we have 

Y H^{W^'\t),t)\j^{t-l), ij{]R\t) 
= Y H^{W^^{t-l),t-l) 

i=l 
ip{M'^,t-l) 

= Y Y,^t_,Ho.iW^'\t-l),t-l), 

i=l 
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the last identity following from the fact that many particles at time t come 
from the same ancestor at time t — 1, with Y^^t-i denoting the number of 
offspring from the i-th particle at time t — 1. 
Integrating on both sides gives that 

^ Ha{W^'\t),t)\T{t-l) 

= Y: E(y)i/,(iy«(i-i),i-i) 

i=l 

i){R'^,t-i) 

i=l 

proving that t i-^ -^Va{t) is a martingale. 

We now show that (^VQ,(t), t > 0) converges to a finite limit almost 
surely. With the above notation we first write 

(2.18) V^{t)= Y: H^(W^^^(t),t)= E Y Ha{W^''"'\t),t) 

1=1 1=1 m=l 

where W^''''^\t) is the m-th child of the Z-th particle which dies at time t — 1. 
Then we can write 

(2.19) = E E Hl{W^^'-\t),t) 

1=1 m=l 
■,p{M'',t-l) Yi,t-i 

1=1 m^n,m,n=l 

+ E E E H^iw^''-^\t),t)H^{w^'-\t),t) 

27^i>i,i=i »"=i n=i 

Therefore 



(2.20) E(Va{tf\j^{t-l), ^{IR\t)) 



tp{M'',t-l) 

i=l 

+ E iy't-i-nt-i)Hiiw(Ht-i),t-i) 

i=l 
ip{M'^,t-l) 

+ Yl Y,,t_^Yj,t_,H^{W^'^{t -l),t- l)H^iW^^\t -l),t- 1). 



Thus 

(2.21) E(Va{tf\j^{t-l) 

= mE{Hl{W^^\t),t)\j^it-l 

i=l 

i=l 
'4){1R'^ ,t-l) 

+ rn^H^{W^\t-l),t-l)H^{W^^\t-l),t-l). 
Y '\mE{Hl{W^'\t),t)\nt-l)] 



i=l 



+ {a^ - m)Hl{W^'\t -l),t- 1)1 + mV,(t - 1)' 



Recall that E(?/'(iR'^,t - 1)) = m^-^ (Athreya and Ney P, p. 9). It is 
easy to see using ^TT^ that E (Hl{W'^^\t),t)) = Hence 



(2.22) E{v^itY) 

= m*-^a!(mtl°l + (a^ - m){t - l)l"l) + m^E (v^it - 1)' 
= m*-^a!(m(tl"l - (t - l)l"l) + ^'(t - 1)1"') + m^E (v^,(t - 1)' 

This gives us that 



(2.23) 



< E 



m 



2t 



m 



2(t-l) 



Hence, using the fact that V^(t)/m* is a martingale we have that 



i=l 



e(E 

(2.24) 



V^{t) V^{t-l) 



m 



t-i 



< E E 



Va(t) Va(t-l)^ 



1/2 



i=l 



E E 

t=i { 

oo 



2\ ^ 1/2 



2t 



m 



2(t-l) 



t|a|/2 

— ^E — ~ ^ 



t=i 



m 



t/2 



SO that 
(2.25) 



E 

t=i 



V^it) V^{t - 1) 



< oo, a.s. 



This shows that (^ 14(i), ^ > 0) converges to a finite limit almost surely. <C> 
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Remark 2.7 Note that by induction from ()2.22p 

t-i t 

(2.26) E(Va{ty) = m*-^a!((T2 ^ m*-^j l°l +m^m*-^(j l°l - (j - 1)1"!)) 

j=i j=i 

and therefore 

oo oo 

(2.27) E (at^) = m-^a\{a^ "^"^i'"' + mJ2 m-^ (jl"l - (j - !)'"!)). 

Lemma 2.8 Let a G 2'^, anc? let Va, Na be as in Lemma \2. (A Then for any 
e > Q, we have that almost surely as t ^ oo, 

(2.28) Ml = iV, + o(m-*/(2+^)). 
Proof. We claim that 

(2.29) ^ = E f ^ m] + o , a.s. 

To see this, we first observe that by Fact 12.41 Chebyshev's inequahty and 
the Borel-Cantelh lemma that almost surely for t ^ cxd, 



max 



X{y,e)-Y.i\{y,e)\T{t) 



Assembling this estimate with ()2.1(jp and Lemma f2.5^ together with the fact 
that sup||y||<i2{i+.) H^iy.t'^) < ct2{i+=)l°l, we get ^HTM . 

Since E(^5^ | r{t)) = ^ (by Lemma it follows from ^TI^ that 



^ -Ml=o fm-*/(2+2^) 1 , a.s. 



As a consequence, 



(2.30) A.„-M)=£^^""^""' ^'^'"y 



This proves our lemma, since e > is arbitrary. <0 
We have now all the ingredients to prove Theorem 11.11 
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Proof of Theorem\r^ Fix A; > 1. Fix < 7 < ^i^, and let t = [T^J. Let 

e > be such that (1 + £:)7 < 2{k+i) ■ ^^^^ show that, almost surely for 
T — > 00, 

(2.31) +0 (r-(^+'^/2)) + O (m-*/(2+.)) ^ 

where V2a-/3 is defined in fl2.16|l . Our Theorem will then follow from Lemma 

EHl 

By Fact 12. 3( we have, almost surely for T ^ 00, 



^^ = ^/ EP^Uy-x)Xiy,t)dx + o{m-'/^'^^^ 
nrJAy ^ 

On the other hand we can write 

(2.32) (27rT)'^/V^^l,(y-x) = ^^J^^y,, exp (-|^) = 
where z = 1/T and 

(2.33) /(., t, ,) = ^^-^ exp (-|^) 

is a C°° function of z near z = as long as ^2; -C 1. If we expand /(z, t, x, y) 
in a finite Taylor series in z around z = 0, it is clear that we can bound the 
remainder Rk+i{z, t, x, y) of order + 1 by a polynomial in \\y — x\\ of order 
at most 2{k + 1). 

According to Lemma 12.51 almost surely for all large T, X{y,t) = as 
long as 1 1?/ 1 1 > T^^^'^^'^ . Together with p.l|) which implies that the number of 
points y with \{y,t) 7^ is bounded by cm* and the fact that A is bounded 
we have 

(2.34) —JY. Rk+i{T-\ t, X, y)X{y, t)dx < cT^^'^^^^^^+^) = o{T). 
rrr J a y 

By inspection of Lemma 12. 2t the first k terms in the Taylor series for 
f{z,t,x,y) give rise to the the first line of ()2.31|1 . completing the proof of 
that formula and hence of our Theorem. ■(> 
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3 convergence 

In this section we show that if the offspring distribution Y has p moments 
for some even integer p then ^^^^ converges in L'p. 
Introduce the notation 

n " 

^n,...,ij=i 5Z 

for summation over non-repeated indices. Let Zt = '^(iR'^,t), Fa-i(t) = 
and 

(3.1) t^aa),...,ab)W=E,t ,=in^aW;..W- 

' ' h=l 

The following Lemma will play an important role in showing that ^^^^ 
converges in L^. 

Lemma 3.1 Let k be an integer with E{\Y\'') < oo. Then for any a^^\ ... ^ a^^^ 
we can find c, (3 < oo independent of t such that 



(3.2) |E(f/,(i)„„,„(.)(t))|<ct^ 



Proof of Lemma \3.1\ We will prove this Lemma by induction on k. The case 
of = 1 is trivial. Assume that we have proven this Lemma for all A; < p — 1. 
We can write 

(3.3) U^w....Mi)=Y.t ,-^\{Fc.i-)■,^,it) 

h=l 

P — Zt-i ^ ( — "n^,. 

fc=i yliu---uylfc=[i,p] h=i \ meAh ^ 

where the sum J2Aiu---uAk=[i,p] ^^^^ over all partitions of [l,p] = {1, . . . ,p} by 
k non-empty sets Ai,...,Ak and Fa-i^mif) = Ha{W^'-'"'\t),t). Introducing 
the falling factorial notation {x)k = x{x — 1) • ■ ■ (x — A; -|- 1)) we have that 



\h=l \ mt^Ah 
k 

nE((i^)^j n ^«('");..(^-i)- 



h=l mSAft 
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Hence 

(3.5) E(f/„,i)^..._,„)(t)|^(t-l 

= E E.tl., E n E n ^.(™);.(^ - 1) 

k=l A-iU---UAk = [l,p]h=l meAh 

= mPf/a(i),...,„(p)(t - 1) 

+EEt:l^, E nE((F)|^j n 

fc=l AiU---UAk = [l,p] h=l m&Ah 

Note that by ^IT[\ 



oo 



rx-tr'^ 12 ^sx-ts^ 12 



(3.6) 5: -ri/n(x, t) E —.Bm[x, t) = e 

— ^{r+s)x-t{r+sf /2 ^trs 

j=0 J- k=0 

Equating coefficients of r^s"^ we find that 

mAn 

(3.7) Hn{x,t)Hm{x,t) =n\m\ V — — p-:Hn+m-2k{x,t). 

^ kl [n — ky\m — k)\ 

Using this to reduce products of Hermite functions to sums we find that 



(3.8) E(f/,a)_„,„(t)|^(t-l)j =m^'f/„,i)^...^,(,)(t-l) 

+ E E c(«; • • • ' - 1) 

i=l/3(i),...,/3(j) 

where Z]/3(i),...,/30) is a finite sum over (3^^\ . . . such that ELi l/^^'^l < 
X]f=i la*-'-*! and the c(a; p; I3^^\ . . . , t) are polynomials in t. Hence by our 
induction hypothesis 

(3.9) E (f/,(i),..vrtW) = "^''E (f/„(i),...,,(p)(t - 1)) +7^„(l),...,„(p)(t) 

with |7^Q,(i),.,,,a(p)(^)l ^ ct'^m^'P^^^^^^^^ for some /?, c < oo independent of t. 
Iterating this completes the proof of our Lemma for k = p. (} 
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Proposition 3.2 Let p be an even integer with E{\Y\^) < oo. Then ^-^^ 
converges in . 

Proof of Proposition W^ Note that because of the presence of the polynomial 



IS 



factor in ()3.2|1 we cannot simply use Lemma 13.11 to show that 
bounded uniformly in LP. Rather, we will show that for some c, /3 < oo 
independent of t 



(3.10) 

Then 
(3.11) 



|E({V;(t) -mK.(t- l)f )| < ct^m*(P-i). 



E 



K(t) v^{t-i)Y^ 



t-i 



< ct^m * 



and therefore (it is here that we need p even) 
V^{t) V^{t-1 



(3.12) 



E 

t=i 



m 



t-i 



< oo 



t=l 



which will complete the proof of the proposition. 

The basic idea of the proof of ()3.1()|1 is that the subtraction eliminates 
the highest order term in the expectation leaving only sums of terms of the 

form f/a(i),...,a('=)(^) ■with k < p — 1. 
We now prove ()3.10|) . We have that 



(3.13) 



E({K(t)-mV;(t-l)f) 



E 

fc=0 

p 

E 



'p' 

V 

, k I 



k=0 

By (jTTH|) we have 



-l)^m'=E {v^-'{t)V^{t - 1)] 
-l)''m''E (e [VP-''{t) I J^(t - 1)) V^{t - i; 



(3.14) 



Zt-i Yi^t-i 

Vait) = E E Fa;l,mit) 
1=1 m=l 



where F^.i^^{t) = H^{W^^^'^\t),t). Thus 
(3.15) V^it) 



E Ej^,...,j^.=i E 

j=l h~\ Hj=n 



n 



n E F'^-.^.At) 



1 "J/ h=\ \ r=l 
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n E f«....w 



h=l \ r=l 



n-1. Zt_i 



n 

, /]_, . . . ,lj 



ih 



iiH h/,=n \ J^' • • • ' "J/ /i=l \ r=l 



and 
(3.16) 



r=l 



I 



h 



s=l 



Thus 



qi + -+qs=lh V^l' ■ • ■ '^V /=1 



(3.17) E ^ F^,,,,.{t)\ :F{t-l) 



r=l 



Ih 



EE((y)j 

91 H hgs=';i 



s = l 



J 



nE(i^.\(t)|^(t-i) 



Using ()3.7|) to reduce products of Hermite functions to sums we find that 
by ()3.15p - ()3.17p we can write, with a^^^ = a, i = 1, . . . ,n 

(3.18) E (C(t) I - 1)) = m"[/„a),...,„(„) (t - 1) 

n-l 

+ J2 c(a; n; P^^^; t)Upw_po) {t - 1) 

where J2/3W,...,pU) is a finite sum and the c{a; n; . . ,p^^^;t) are polyno- 
mials in t. 

We next observe that 

E Ej^^...^j^.=i E 

j=l h-{ \-lj=n 

= ^Q(i),...,a(")(i - 1) 



n 



Ji, . . . ,1 



J/ h=l 



n-l Zt-i 

E Ej, j.=i E 



n 



(iH \-lj=n 

14 



— ^a(i),...,Q(")(i — 1) 
n-1 

+ E E d{a; n- t)U^w_^u) {t - 1) 

i=i 7(1),. ..,7U) 

where we have again used ()3.7p to reduce products of Hermite functions to 
sums, and the d{a; n; 7'^-'^^ . . . , 7*^-'^; t) are polynomials in t. 
Similarly 

(3.20) Upw_f,o) (t - l)U^w_^i,) (t - 1) 

E';:;,.,nv);j^-i)) (El;:i.,n^.<o;,(t-i)) 

= ^/3(1),...,/3(J),7(1),...,7(*)(^ ~~ 1) 

+ E E f{P,rX^'\..-X^'^^;t)u^w,...,c^m,{t-i) 

m=l 

where we have abreviated P = {P^^\ . . . , P^^''), 7 = (7'-"^^ • • • , 7*-'^-'). 
Combining ()3.18|) - ()3.20|1 we have that for each k < p 

(3.21) m'^E (V^'it) I ^(t - 1)) V^{t - 1) 
= mFU^w^_^^(v){t - 1) 

+ E E 7^'\ . • . , 7^^'^; t)f/7W,...,7U) (t - 1) 

i=l 7(1),...,7(J) 

where the /(a; n; 7*^^^ . . . , 7*^-'^; t) are polynomials in t. Substituting back 
into ()3.13|) and using the fact that YX=q = we find that the 

m*'[/Q,(i)^ „(p)(t — l)'s cancel, and we can write 

(3.22) E({K(t)-mK(t-l)f) 

= E E ^?(«; ^; 7^'\ • • • , 7^^'^; m {u,w,...,yo^ (t - 1)) 

i=l 7(l),...,7(i) 

where the (7 (a; n; . . . , 7^-'^; t) are polynomials in t. ()3.2|1 then completes 
the proof of ()3.10|) and hence of our Proposition. <(> 

Remark 3.3 Note that by Proposition 13. 21 we have that ||-^^^^||p is bounded 
uniformly in t, so that 

(3.23) \Km<om'. 
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Arguing as before, any f/Q(i),...,Q,(p)(t), where a^^\ . . . , a'^'^^ are now arbitrary, 
can be written as 



p 

E 

i=l 



(3.24) U^(i),...,aip) (t) = n V^i.) (t) + terms of 'lower order' 



and thus using ()3.23|1 . Holder's inequality and ()3.2|) for < p — 1 we can 
refine ()3.2j) and find c, /3 < oo independent of t such that 



(3.25) |E(t/„(i),...,,,.,(t) 



(Here we require that Y have r momnets for some even r > p). 
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